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Thermal Analysis on Cryogenic Liquid Hydrogen 
Tank on an Unmanned Aerial Vehicle System 
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Thermal analyses are performed on the liquid hydrogen (LH2) tank designed for an un- 
manned aerial vehicle (UAV) powered by solar arrays and a regenerative proton-exchange 
membrane (PEM) fuel cell. A 14 -day cruise mission at a 65,000 ft altitude [ 1 ] is consid- 
ered. Thermal analysis provides the thermal loads on the tank system and the boiling-off 
rates of LH2. Different approaches are being considered to minimize the boiling-off rates 
of the LH2. It includes an evacuated multilayer insulation (MLI) versus aerogel insulation 
on the LH2 tank and aluminum versus stainless steel spacer rings between the inner and 
outer tank. The resulting boiling-off rates of LH2 provided by the one-dimensional model 
and three-dimensional finite element analysis (FEA) on the tank system are presented and 
compared to validate the results of the three-dimensional FEA results. It concludes that 
heat flux through penetrations by conduction is as significant as that through insulation 
around the tank. The tank system with MLI insulation and stainless steel spacer rings 
results in the lowest boiling-off rate of LH2. 


I. Introduction 

The densified propellant (LH 2 ) is stored at cryogenic conditions near its boiling point. Any heat trans- 
ferring into the propellant could cause phase change such that the storage of the propellant becomes more 
unmanageable and unpredictable. Heat leaks into the storage tank must be minimized through proper in- 
sulation. Boiling-off rate of the propellant in the storage system has to be assessed. For a 14-day mission, 
minimizing the boiling-off rate of LH 2 could reduce the system mass significantly. Both an evacuated MLI 
blanket and aerogel insulation are investigated for comparison. An evacuated MLI blanket gives better 
thermal insulation; however, the MLI fabrication and installation requires extreme care and precision. In 
addition, the thermal performance of MLI is very sensitive to the vacuum pressure [2]. 

A representative LH 2 tank system is shown in Fig. 1. It consists of two concentric spherical tanks 
connected by spacer rings, which serve as a sleeve bearing for the support rod passing through the tank 
assembly. The inner tank is welded to the support rod, and the two rings are welded to both inner and outer 
tanks to seal the vacuum gap between the tanks where either MLI or aerogel insulation is used. The inner 
tank has an 8.5 ft diameter with a 0.031 in. thickness and holds 1,300 lbm of LH 2 stored at an internal 
pressure of 30 psia with 5 % ullage. The outer tank is 0.25 in. thick and the MLI insulation is approximately 
2 in. thick. The support rod has a 4 in. inner diameter with a 0.1 in. thickness while the ring has a 4.5 in. 
inner diameter with the same thickness. The support and inner/outer tanks are made of 2014— T6 aluminum 
and the spacer ring material is either 2014— T6 aluminum or stainless steel for comparison of heat loss. 

In the following, the one-dimensional thermal circuit for the tank system is presented first, which is 
followed by the three-dimensional FEA on the tank system and conclusions. 
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(a) Cutaway view (b) Tank/rod interface detail 


Figure 1. A representative LH 2 tank design (courtesy of Timothy Roach). 


II. One-dimensional Thermal Circuit for the Tank System 

The heat transfer modes across the tank wall are sketched in Fig. 2, showing free convection and radiation 
outside the outer tank and free convection inside the inner tank. Note that the wall thickness of both the 
inner and outer tanks is ignored due to the relatively thin wall. Between the inner tank and outer tank, 
heat flux transfers through three parallel paths, which shows conduction through two rings and MLI/aerogel 
insulation. Across the evacuated MLI or aerogel insulation, both radiation and conduction take place, which 
is not shown in Fig. 2. Either the effective thermal conductivity, k e , or effective heat transfer coefficient, 
/i e , can be used to represent the evacuated MLI or aerogel insulation [3]. The equivalent thermal circuit is 
plotted in Fig. 3. T a and Tlh 2 are the temperatures of the ambient air and liquid hydrogen, respectively. 
The temperature of the inner and outer tanks is T) and T 0 , respectively. Assuming the radius of the inner 
and outer tanks being ri and r 2 , respectively, the thermal resistance between the outer tank and ambient, 
R\ ((W/K) -1 ), is defined as 

Ri = [1/(470-2*!) + l/(47rr 2 2 ea(T 0 + T a )(T 0 2 + T 2 ))] (1) 

where h\ being the heat transfer coefficient (for convection mode), e being the surface emissivity(for radiation 
mode), and a being the Stefan-Boltzmann constant (a = 5.67 x 10 -8 W/m 2 -K 4 ); the thermal resistance 
across the insulation annulus (MLI/aerogel), i ?2 ((W/K) -1 ), is defined as 

R -2 = l/(4nk e )(l/r 1 - l/r 2 ) (2) 
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T a = 216.7 K at 65,000 ft altitude 
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Figure 2. Schematic view of the one-dimensional heat transfer across the tank wall. 


r 3 



the thermal resistance through the conduction (outer tank, ring, inner tank, and support), R 3 ((W/K) x ), 
is defined as 

B 1 | 1 | 1 | 1 

k\A c \jL\ kiAcpjLi k^Aa^/L-s k^Ac^/L^ 

where k n ,L n ,A c n with n = 1,2,3, and 4 are the thermal conductivity, the length of conduction, and the 
contact surface area for the outer tank, the ring, the inner tank, and the support rod, respectively. The 
thermal resistance between hydrogen and the inner tank, R 4 ((W/K) -1 ), is defined as 

Ra = 1/(47 nihi) (4) 

with /14 being the heat transfer coefficient. The total thermal resistance of the circuit, R ((W/K) -1 ), is 

R — Ri + Ra 4 — j 2 ~ (5) 

lh + R ~ 3 

The heat flux, q (W), can be computed by using 

<7 = (T a ~ T LH2 )/R (6) 

and the temperature of the inner and outer tank are 

T{ = Tlh2 + qRa, T 0 =T a - qRi (7) 
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The heat flux transferring through R -2 and Rj , denoted as q> and q$, respectively, are computed as 


92 — (T 0 — Ti)/R‘i. 93 — (T 0 — T)/i?3 
The boiling-off rate of LH 2 , M (kg/s), is computed by using 


(8) 


M = q/hf g 


( 9 ) 


Note that the heat flux to LH 2 prior to phase change is ignored since the amount of heat is very small 
compared to the latent heat. 

At a 65000 ft altitude, the ambient air temperature is —70° F (216.7 K) and the pressure is 0.826 psia. 
The LH 2 is filled in at a temperature of 20.4 °K. At a 30 psia internal pressure, the boiling temperature of 
LH 2 is 23 °K and the latent heat, hf g , is 428 kJ/kg. In Eq. (1), hi is chosen as 1 W/m 2 -K by using the 
equation for free convection around a sphere [4], and e = 0.02 for polished aluminum surface. In Eq. (2), 
k e = 0.00016 W/m-K for MLI and 0.0017 W/m-K for aerogel [4]. In Eq. (3), = 12.6 W/m-K for stainless 

steel AISI 304 and 163 W/m-K for aluminum, L 2 = 2 in. and A c> 2 = 0.44n in. 2 , and Aq = k^ = Aq = 163 
W/m-K for aluminum inner and outer tanks and support rod. In Eq. (4), At is chosen as 10 W/m 2 -K based 
on the table data for free convection of liquid [4]. 

In the first case, it is assumed that no heat leak through the ring by conduction, that is, R 3 = 00 . The 
heat flux and thermal resistance of the tank system is listed in Table 1. It can be seen that the thermal 
resistance between the hydrogen and inner tank and that between the ambient air and outer tank are very 
small. Replacing MLI with aerogel insulation, the thermal resistance, R 2 , increases almost 10 times. The 
total heat flux is also increased almost 8 times. In the second case, it is assumed that the ring is made of 
stainless steel, the thermal resistance of the ring becomes dominant in R 3 and the resistance of aluminum 
tanks and support is negligible, which gives 


r 3 


1 

feA C] 2/T2 


(10) 


The calculated results of the thermal resistances and heat flux for both MLI and aerogel insulation are listed 
in Table 2. It can be seen that the total heat flux is increased to 97.8 W from 13.3 W for MLI and to 217 W 
from 111 W for aerogel insulation. It is obvious that the heat loss through the penetration by conduction is 
significant. The total heat flux increases 7.35 times and 86% heat loss through the rings for MLI and 2 times 
and 38% heat leaks through the rings for aerogel. For the tank system with aluminum rings, the ID result 
is not available due to the thermal resistances of the inner/outer tanks and supports can not be estimated, 
which will be investigated in the three-dimensional thermal analysis as follows. 


III. Three-dimensional FEA of the Tank System 

The tank system with MLI or aerogel insulation and aluminum or stainless steel spacer rings are analyzed 
using MSC Pthermal. A part of the finite element model of the tank system that has 73,494 nodes and 213,955 
tetra elements is shown in Fig. 4. The free convection boundary condition is imposed on the inner surface 
of inner tank and the outer surface of the support inside the tank. The heat transfer coefficient is chosen as 
5 W/m 2 -K for the gaseous hydrogen part and 10 W/m 2 -K for the liquid hydrogen part based on the table 
data for free convection in [4]. The part of the support outside the tank is assumed to be adiabatic. Free 
convection and radiation boundary conditions are imposed on the outer surface of the outer tank in the same 
way as that in one-dimensional analysis. On the outer surface of the inner tank, the effective heat transfer 
coefficient, h e = k e /t with t being the thickness of the insulation, was used, and h e = 0.00314 W/m 2 -K for 
MLI and 0.03346 W/m 2 -K for aerogel. There is no contact between the support and outer tank. In the case 
with the rings and support rod, the support and inner tank, the rings and inner/outer tanks are assumed to 
have no thermal resistance between each of them. Otherwise, larger thermal resistances are imposed at the 
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Figure 4. Part of the finite element model of the tank system. 


interface to represent no heat transfer through there. Three cases including (i) without rings and support; 
(ii) with support and aluminum rings; and (iii) with support and stainless steel rings are considered here for 
both MLI and aerogel insulation. 

To show the representative temperature distribution on the tank assembly, the MSC Pthermal results 
of the tank system with evacuated MLI blanket and aluminum spacer rings are plotted in Figs. 5 to 12. 
The temperature contour for the tank system is shown in Fig. 5, showing that most part of the outer tank 
is almost as hot as the ambient and most of the inner tank is nearly as cold as the hydrogen. A large 
temperature gradient appears around the top and bottom areas where the support penetrates. The close 
view of the top area is shown in Fig. 6. It can be seen that the heat conducts from the hot outer tank 
to the connecting ring and the inner tank then goes to the support. The hottest spot on the top part of 
the support is located at the interface of the support and inner tank, showing a temperature of 126 °K. 
The similar phenomena is seen at the bottom part plotted in Fig. 7. However, the bottom part has a 
lower temperature due to the contact with liquid hydrogen instead of gaseous hydrogen. The maximum 
temperature on the bottom part of the support is around 114 °K. The top ring is around 10 °K hotter 
than the bottom ring. The observation described above is further verified by the plots of the temperature 
contour on the top/bottom ring (together with part of the inner and outer tanks) in Figs. 8 and 9, and on 
the support (together with part of the inner tank) in Fig. 10. The temperature contour of the inner and 
outer tanks are plotted in Figs. 11 and 12, respectively, to show a clear view of the temperature distribution. 
The corresponding temperature contours for the tank system with aerogel insulation are plotted in Figs. 
13 to 15. The temperature distribution is very similar to those presented above for the tank system with 
evacuated MLI. However, the outer tank is much hotter and the inner tank is much cooler because much 
less heat transfers from the outer tank to the inner tank. The temperature distribution of the tank system 
with stainless steel spacer rings is very similar to that with aluminum spacer rings and is not plotted here. 

The three-dimensional results of the total heat flux are listed and compared with the one-dimensional 
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result in Table 3 for the three cases. For the system without a support and rings, the heat flux transferring 
into the LH 2 is 16.5 W for MLI and 137 W for aerogel, which agrees well with the one-dimensional results 
that are 13.3 W and 111 W, respectively. For the system with stainless steel rings, the total heat loss is 87.5 
W for MLI and 213 W for aerogel, which is very close to the one-dimensional results that are 97.8 W and 
217 W. For the system with aluminum rings, the total heat flux (237 W) for MLI is 2.7 times of that with 
stainless steel rings, and while for aerogel, the heat flux (354 W) becomes 1.7 times. It can be concluded that 
the system with stainless steel spacer rings and MLI blanket results in the least heat flux. The corresponding 
boiling-off rate of LH 2 , M, is 1.62 lbm/hr. For a 14-day mission, an extra 544 lbm mass of LH 2 has to be 
carried for boiling-off. However, the tank is designed to carry 1,300 lbm LH 2 for the required power. With 
the total boiling-off mass of the LH 2 within the 14-day mission being taken into an account, the vehicle can 
operate 8 to 9 days on 1,300 lbm LH 2 instead of 14 days. The heat leak from the tank system with stainless 
steel rings and MLI blanket needs to be further reduced to minimize the boiling-off rate of LH 2 . A 81.1% of 
total heat loss still goes through the penetration by conduction. The ring between the inner and outer tank 
could be replaced with a bellows made of stainless steel to increase its thermal resistance. The support rod 
could be modified to have two separate parts located at the top and bottom parts of the tank and ended at 
the inner tank such that the heat transfer area between the support and hydrogen is greatly reduced. Other 
solutions that might be pursued to minimize the heat flux transferring to the hydrogen from the ambient 
are not the focus of this paper and will not be further discussed here. 

IV. Conclusion 

The thermal analyses on the LH 2 tank system used in an unmanned aerial vehicle were performed. Both 
one-dimensional thermal model and three-dimensional FEA results are presented and discussed. The com- 
puted results of the heat flux transferring to the tank system showed that the heat leak through penetrations 
are significant. For a 14-day mission, the current tank design with MLI insulation requires an extra 544 
lbm LH 2 for boiling-off while the tank carries only 1,300 lbm LH 2 for the required power. The vehicle can 
operate 8 to 9 days when taking into an account of the boiling-off of the LH 2 . The heat flux transferring to 
the tank system from the ambient has to be further reduced. The thermal resistance between the inner and 
outer tanks should be maximized by reducing the contact surface area between the ring and the inner/outer 
tank and by replacing the ring with a bellows made of stainless steel if it is possible. 
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Table 1. One-dimensional thermal results (without support and rings) 
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Table 2. One-dimensional thermal results (with support and stainless steel rings) 
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Table 3. Three-dimensional thermal results compared with the one-dimensional results 
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Figure 5. Temperature (K) contour on the tank system with MLI. 
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Figure 6. Close view of the top part of the tank system with MLI. 
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Figure 7. Close view of the bottom part of the tank system with MLI. 
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Figure 8. Temperature (K) contour of the top ring with part of the inner /outer tank (MLI). 
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Figure 9. Temperature (K) contour of the bottom ring with part of the inner /outer tank (MLI). 
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Figure 10. Temperature (K) contour of the support rod with part of the inner tank (MLI). 
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Figure 11. Temperature (K) contour of the inner tank (MLI). 
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Figure 12. Temperature (K) contour of the outer tank (MLI). 
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Figure 13. Temperature (K) contour on the tank system with aerogel insulation. 
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Figure 14. Temperature (K) contour of the inner tank (aerogel insulation). 
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Figure 15. Temperature (K) contour of the outer tank (aerogel insulation). 
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